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Abstract 

We refer to the ground state of a gravitating, charged ideal fluid of fermions held at 
a finite chemical potential as an 'electron star'. In a holographic setting, electron stars 
are candidate gravity duals for strongly interacting finite fermion density systems. Wc 
show how electron stars develop an emergent Lifshitz scaling at low energies. This IR 
scaling region is a consequence of the two way interaction between emergent quantum 
critical bosonic modes and the finite density of fermions. By integrating from the IR 
region to an asymptotically AdS4 spacetime, we compute basic properties of the electron 
stars, including their electrical conductivity. We emphasize the challenge of connecting 
UV and IR physics in strongly interacting finite density systems. 



1 The broader context 



A challenge facing contemporary condensed matter theory is the description of a 2+1 di- 
mensional finite density of fermions interacting with a gapless collective bosonic excitation, 
such as a spin density wave or emergent gauge field. Such theories arise, for instance, when 
a Fermi liquid is tuned across a quantum phase transition. The low energy dynamics of 
the system of fermions interacting with the critical bosonic mode can be characterised as 
metallic quantum criticality. While in 3+1 dimensions one can proceed to integrate out 
the fermions and obtain a stable Gaussian theory for the boson [I], this approach does not 
give correct answers in 2+1 dimensions, see e.g. [11 OH], as it ignores an infinite number of 
nonlocal marginal couplings in the effective theory for the boson. One should not integrate 
out the fermions in this case but rather flow to a scaling regime involving both the boson 
and fermion fields. The resulting low energy theory is strongly interacting, e.g. [3]. 

One might have hoped to perform a (vector) large N analysis as a perturbative handle 
on the theory. It has recently been demonstrated El H] that the vector large ./V expansion 
breaks down for 2+1 dimensional metallic quantum critical systems. This occurs because a 
potential IR divergence at high loop order is cured by a self-energy of order 1/N, leading to 
extra factors of N in the numerator in certain Feynman graphs. Partially motivated by these 
difficulties, in this paper we will use the holographic correspondence [TJ 13 El QUI to study 
a strongly interacting system of gapless bosons with a finite density of fermions. Before 
proceeding we should note that more traditional approaches to this problem have also been 
proposed |11| H2l [13] and that our framework does not appear to include ingredients that 
are likely crucial for applications to the original systems of interest, such as Fermi lines with 
cold regions as well as hot spots. We will, however, describe the emergence of a low energy 
scaling regime from the interaction of critical bosons with a finite density of fermions. The 
essential physics of this process was noted in [14] • 

In the holographic correspondence a charge density is implemented by a bulk Maxwell 
field, dual to the current operator in field theory. The asymptotic boundary value of the 
Maxwell field determines the chemical potential of the field theory. This is a UV input 
the consequences of which we wish to explore at low energies. In the simplest bulk setup 
of Einstein-Maxwell theory, the gravitational solution dual to the finite chemical potential 
theory is then determined uniquely to be the planar AdS-Reissner-Nordstrom black hole. 
This black hole was therefore a natural starting point for investigations into strongly inter- 
acting finite density systems [151 EH El ES] ■ Reviews of this and other earlier work can be 
found in $ E1H3I. 



1 



A conceptually problematic aspect of charged black holes in an applied holography 
context is their blackness. One is often interested in temperatures much lower than the 
scale set by the charge density. In this extremal black hole limit the horizon remains 
present and the actual source of the bulk electrical field remains hidden. Thus, within 
a bulk effective field theory approach to holography we don't have explicit access to the 
zero temperature charged degrees of freedom^] Although the universality of a black hole 
description of charge density is appealing, and may yet have important consequences, it 
makes it difficult to connect with basic experimental implications of a finite density that 
depend on the nature of the charge carriers, such as a Fermi surface in the case of fermions. 
A Fermi surface appears not to be inherent to the gravitational geometry, but depends on 
the nature of external probes [2T1 1221 12511251 125] . 

It is perhaps fortunate, therefore, that low temperature charged AdS black holes are 
found to be unstable towards a range of processes that discharge the black hole and can 
lead to spacetimes without black hole horizons. The instabilities include condensation of 
charged scalar fields [26J, Cooper pairing of charged fermions [27], emission of D branes 
[28\ [2U1 [T3] , backreaction of a bulk fermionic charge density induced by the local chemical 
potential [H] , confinement [30l |3_TJ [32] , and perhaps the emergence of underlying lattice 
degrees of freedom |33j . It is not clear at this stage whether all zero temperature charged 
AdS black holes with a finite size horizon are unstable |34j . If they are, this fact may be 
closely tied up with a version of the 'weak gravity' conjecture |35j. The instabilities lead to 
a new zero temperature bulk geometry, often without a finite size horizon (e.g. [36, 37J). 
The charge is then carried by explicit bulk fields and we can identify the corresponding field 
theory operators as responsible for the finite density dynamics. 

A fruitful approach taken in previous works in order to explicitly model holographic 
charge carriers even in the presence of horizons is to add probe D branes into the bulk, see 
e.g. |38[ l39l l4"0l 114] . The limitation of this approach, shared with that of probe fermions 
[21~1 122| [23l [24] , is that it does not capture the two way interaction between the (putatively 
fermionic) charge carriers and the quantum critical modes. 

In this paper we will expand upon section 7.4 of [14J and describe the electromagnetic 

and gravitational backreaction of charged fermions on the holographic spacetime geometry. 

In general this is a very difficult problem as the fermions cannot be treated classically. The 

1 While in some supersymmetric theories we might hope to be able to adiabatically continue the problem 
to a weakly coupled regime and 'count' the degrees of freedom there [2D], this does not help us with our 
objective of understanding the strongly interacting finite density dynamics on its own terms. 
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coupled fermion-Maxwell-gravity system becomes tractable in a limit in which the fermions 
may be treated locally in the bulk as an ideal fluid of zero temperature charged free fermions. 
This approach mirrors the standard description of neutron starts in astrophysics, following 
the original Oppenheimer-Volkoff-Tolman papers |414 [22] . The neutron star equations were 
generalised to an asymptotically AdS setting in [33]. Given that our fermion fluid is charged, 
we will refer to our solutions as electron stars. 

In the following section we set up the equations of motion for a charged ideal fluid in 
Einstein-Maxwell theory with a negative cosmological constant. Later in section [4] we derive 
these equations from an action. We discuss the regime of validity of the fluid description 
depending on the Newton and Maxwell couplings as well as the cosmological constant. In 
section [3] we characterise the (planar) electron star solutions to these equations of motion. 
We show that the IR of the geometry has an emergent Lifshitz scaling and compute the 
dynamical scaling exponent z as a function of the parameters of the theory. Integrating 
out from the IR region to the spacetime boundary, we numerically obtain the full electron 
star solutions and compute their mass and charge. By perturbing the solutions in section [5] 
we obtain the electrical conductivity as a function of frequency. We find that the electron 
star conductivity exhibits a universal low frequency behaviour previously noted in other 
solutions with an IR Lifshitz scaling. Our work presents a framework in which the physics 
of a strongly interacting finite fermion density system can be investigated; we enumerate 
some of the more pressing open directions in the final discussion section. 



2 Equations of motion: background 

We are interested in 3+1 dimensional zero temperature configurations of a charged perfect 
fluid in a holographic setting. We will introduce an action principle in a later section, but 
for the moment will work with equations of motion. The Einstein-Maxwell equations with 
a negative cosmological constant and sources are 

Rab - \g ab R - ^9ab = k 2 (F ac F b c - - A g ab F cd F cd ^ + T ab \ , (2.1) 

and 

V a F ba = e 2 J b . (2.2) 
Here the perfect fluid energy momentum tensor and current are 

Tab = {p + p)u a u b + pg ab , J a = au a . (2.3) 



3 



The four velocity u should be normalised so that u 2 = — 1. The second Bianchi identity 



requires that the right hand side of (2.1) be transverse. Similarly (2.2) requires that the 
current be conserved. In the above expressions the cosmological constant scale L, Maxwell 
coupling e and Newton constant k 2 are constants while the pressure p, energy density p and 
charge density a are fields on spacetime that will be related through the equation of state 
of the fluid. 

For the background we wish to make the following 'planar star' ansatz for the metric 
and Maxwell field 

ds 2 = L 2 (-fdt 2 + gdr 2 + ^ (dx 2 + dy 2 )^j , A=^hdt. (2.4) 

Here f,g,h are functions of the radial coordinate r. The pressure and energy and charge 
densities are also functions of r. It is useful to scale out the couplings and write 

p= zk? p > P= Z^' a = ^?- (2 - 5) 

The velocity has nonzero component u l = \/{L\ff). 

It is straightforward to show that the above Einstein-Maxwell equations are solved 
provided that the following four equations are satisfied 

p' + (p + p)^-^j = 0, (2.6) 
1 1 1 + - + -) +(p + p)g = , (2.7) 



/ 9 r 



+ g(3+p)-^ = 0, (2.8) 



ti' + ~g(p + p)-gVfZ = 0. (2.9) 

These are four equations for six variables and so an additional equation of state must be 
specified in order to close the system. One of the equations is second order. While we 
could set h! = F at this point to obtain purely first order equations in terms of the Maxwell 
field strength, we are shortly about to include the effects of Thomas-Fermi screening which 
introduces an explicit dependence on the Maxwell potential h. 

We will focus in this paper on the case in which the ideal fluid is made from zero 
temperature charged fermions with mass m. Firstly recall that in flat 3+1 dimensional 
space with chemical potential p we would have 

p= I Eg(E)dE, a= I g{E)dE , -p = p-p J a. (2.10) 

Jm J m 
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The last of these expressions is the usual thermodynamic relation for the grand canonical 
ensemble. We have taken the charge of the fermion to be one in units where the Maxwell 
action is \F 2 . The density of states is 



g(E) = f3E^fW 



(2.11) 



The constant of proportionality ft is order one, the exact value is not important for us. We 
will see shortly that it is a rescaled constant /3 that we wish to dial. Finally, if [x < m then 
no states above the vacuum are populated and so p = p = a = 0. 

We will work in the approximation in which the fermion physics is determined by the 
local chemical potential, which is the tangent frame value of the background Maxwell field 



e h 



(2.12) 



V7 

This 'locally flat space' approximation will be shown to be self consistent in an interesting 
regime of parameters shortly. In section [4] we will derive the relation (2.12) from the same 
action that implies the ideal fluid-Einstein-Maxwell equations of motion. Substituting into 



the flat space formulae ( 2.10 ) and scaling the integration variable leads to the 'dimensionless' 
expressions 
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(2.13) 



(2.14) 



Again, the energy density and other variables vanish if -4^ < rh. The integrals in (|2.13|) are 



easily performed analytically. The local free fermion equation of state described by (2.13) 



does not include corrections due to gravitational and electromagnetic interactions. We will 
check below that these corrections are negligible in the regime in which we will work. 



The ansatz (2.13) determines three of our six functions, and therefore we have to check 



that it is consistent with the four equations (2.6) - (2.9). Indeed this is the case. The first 



equation (|2.6|) is in fact closely related to the first law of thermodynamics and is satisfied 

.uce t 
gha 



by (2.13). The four equations of motion then reduce to the following three equations 

7' 
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f 



g' 4 
, +- + - 

f 9 r 
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(2.15) 



In these expressions p and a are given by (2.13). 

Before solving these equations, we should discuss the values of the two free parameters 
(3 and m 2 . Recall that in the classical gravity regime n/L <C 1. We will see shortly that 
the interesting regime we wish to explore in this paper has the 'scaled' constant /? of order 
one. In order to achieve this we therefore need e 2 ~ n/L <C 1. Curiously, this is a fairly 
natural relationship from the point of view of string theory, as it requires the gravitational 
('closed string') coupling to be the square of the Maxwell ('open string') coupling. This 
usually corresponds to the 'probe brane' limit; it is interesting that integrating out fermions 
charged under a probe brane gauge field results in an order one local backreaction on the 
spacetime, in the 'dimensionless' sense that we mean it. 

In the regime in which /3 is order one, the dimensionless mass squared is then of order 
m 2 ~ e 2 m 2 L 2 . The operator dual to the bulk fermion might typically have a scaling 
dimension A ~ mL ~ eL/n ~ 1/e 1, leading to rh 2 ~ 1. For the moment we will 
therefore take rh to be order one (including m = 0). 

At this point we can check whether the approximation of using the local flat space 
results is valid. One requirement for the flat space treatment is that the density of fermions 
is large compared to the curvature scale of the geometry. Thus we can compute, under the 
assumption that a is order one and that e 2 ~ k/L <C 1, 



L 1 

crL 3 ? >1. (2.18 

eK e 4 



Therefore the regime of order one backreaction of the fermions together with the classical 
gravity limit is compatible with our 'Thomas-Fermi' treatment of the fermions. Also com- 
patible with our treatment is the fact that mL 3> 1, implying that the Compton wavelength 
of the fermions is much smaller than the curvature scale (in the massless case one can note 
that fi\oc.L 3> 1). 

We can also now check the validity of our 'mean field' description of gravitational and 
Maxwell interactions. In this description, the interactions are between local charge and 
energy densities, but the equation of state determining these densities does not incorporate 
these interactions. Following [33] we can estimate the local effect of interactions through 
the Boltzmann formula 

^-aVfl. (2.19) 

Here vf is the order one Fermi velocity and f2 is the gravitational or Maxwell total cross 
section. Using the scaling of various quantities given above, we can easily estimate that 
the dimensionless quantity {a ri\ oc )~ l da / dt ~ e 4 < 1 for both gravitational and Maxwell 
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interactions. Thus the local effect of interactions is parametrically negligible. We now 
proceed to solve the equations of motion (2.15) - (2.17) treating (3 and rh as order one free 
parameters. 



3 Solution to the background equations of motion 

3.1 Low energy scaling regime 



In the IR region of the geometry, which will be r — > oo in our coordinates (2.4), one 



finds [14J an emergent Lifshitz scaling [44]. In fact the Lifshitz metric is an exact solution 



to the equations of motion (2.15) - (2.17). This is perhaps intuitively reasonable: the 



effect of having a local charge density given by the local background chemical potential 



(2.12) is to screen the electric field. This might be thought of as a form of Thomas-Fermi 
screening. Once the electric field has a 'mass' it cannot support an AdS2 extremal near 
horizon geometry. Instead, massive vector fields are known to give rise to Lifshitz solutions 
|44j . The metric and Maxwell functions take the form 

/ = ^I> 9 = ^2, h = ~T- 

Here z is called the dynamical critical exponent and is given in terms of j3 and rh by plugging 
the above Lifshitz ansatz into the equations of motion. From two of the equations of motion 
we find 

hl = — , 9l= 36 ^-^ 4 , . (3.2) 
z ' y °° ((l-m 2 )z-l)3/3 2 

The remaining equation of motion then gives a complicated relationship between z, rh and 
P which we cannot solve explicitly. The dependence of z on f3 in plotted in figure [l] below 



for three values of rh. Note that the local chemical potential (2.12) is constant on these 
backgrounds. 

It is possible to extract the asymptotic behaviours analytically. For fixed rh at large 
(3 — > oo: 

1 6 4 / 3 m 2 / 3 1 

Z ~ 1 - rh 2 + (rh 2 - l) 4 /3(2m4 _ 7rh 2 + 6) 2 / 3 /3 2 /3 + " ' ' ( ' ' 

In the massless limit the expansion is a little different 

6 

z = l + - + •••, (f3^oo,rh = 0). (3.4) 
$ 

At fixed rh and small /3 — > (this was the limit considered in 



36 1 3m 4 log i+VT^? 



1 + n , ° 9 " + • • • ■ (3.5) 



(1-^2)3/2^ 2(1 -m 2 ) 3 / 2 
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Figure 1: Dependence of the IR dynamical critical exponent on /3. From left to right, the 
three curves have rh = 0,0.55 and 0.7. 

It is immediately seen that these asymptotic results agree with the behaviour exhibited 
in the above figure. We see that at intermediate values of (5 the dependence of z on (3 
interpolates between the limiting behavious without any intermediate features. 

To make sense of the above expansions, we should first note that if z — > oo the geometry 
becomes AdS2 x M 2 . We see that this occurs as rh — > 1 from below or as (3 — > at fixed 
rh. In these limits the fermion backreaction is being turned off and one recovers the near 
horizon geometry of an extremal planar Reissner-Nordstrom-AdS black hole. In particular 
it is clear that the above solutions only make sense for 

0<m<l. (3.6) 

For masses bigger than unity, the Lifshitz background chemical potential is not able to 
induce a density of fermions. It is possible that interesting scaling behaviour arises in the 



limit rh — > 1, but we will not investigate this here. From (3.3) and the above plot we can 



see that by dialing (3 at fixed rh we can achieve all z satisfying 



* > ^ 1 • ( 3 - 7 ) 

1 — m z 



It is also interesting to note that for massless fermions z — > 1 as /3 — > oo, hence the 
geometry becomes AdS^. In general, the emergence of an IR Lifshitz scaling geometry with 
dynamical critical exponent z tunable using couplings and the mass is reminiscent of similar 
results for holographic superconductors [36l EZ] • The physical difference between the two 
cases is that for the superconductors the bulk Maxwell field becomes massive due to the 



Anderson-Higgs mechanism while in the present case the 'mass' is due to screening by the 
charge density. A different type of screening of the Maxwell field, due to a dilaton coupling 
rather than a charge density, was shown to lead to a Lifshitz IR region in \A5\ 06] . 

The IR Lifshitz solution has the dual field theory interpretation of a low energy scaling 
regime arising from the interaction of a finite density of fermions with emergent critical 
bosonic modes (the metric and Maxwell fields in the IR of the bulk). To explicitly connect 
this scaling to the presence of a finite charge density, we need to integrate out to the UV 
boundary of the spacetime, where the charge density appears as a boundary condition. 
We do this in the following subsections. The presence of such density-induced emergent 
quantum criticality is a non trivial and phenomenologically exciting aspect of our models. 



3.2 From the scaling regime to the electron star boundary 

The next step is to flow up the holographic renormalisation group flow. We do this by 
starting with the Lifshitz IR fixed point of the previous subsection and perturbing it by an 
irrelevant deformation. We then follow the flow induced by this deformation into the UV 



by (numerically) solving the differential equations of motion (2.15) - (2.17). Our treatment 
here is very similar to that of |36} 05] . The main difference with those works is that the 
electron star will 'end' at some specific radius r s where the fluid pressure and charge and 
energy densities all go to zero. 

To perturb away from the scaling solution we can write 

/ = ^(l + /ir a + ---), g= 9 -^(l + gir* + ...), h= h ^(l + h 1 r" + ...). 

(3.8) 

We are looking for solutions where the perturbation grows towards the UV (r — > 0) and dies 
off in the IR (r — > oo). By substituting the above expansion into the equations of motion 
one easily finds that the three allowed exponents are 

2 + z J9z 3 - 21z 2 + 40z - 28 - m 2 z(A - 3z) 2 , . 

a ° = 2 + Z ' a± = — ± 2V(l-m')*-l • (3 - 9) 



As in [361115], the 'universal' relevant deformation with exponent q.q presumably generates 
the finite temperature solution, which we will not consider here. The two exponents a± 



are both real for the range (3.7) of z and m that we have access to. The exponent a_ is 
negative and therefore this is the mode that we need to follow. The q_ and a+ modes 
correspond, respectively, to the coupling and expectation value of an irrelevant operator in 
the IR theory. At a practical level, the presence of the ceo and a+ modes, which must be 
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set to zero for a regular IR Lifshitz region, is why one must numerically integrate from the 
IR outwards towards the boundary rather than the other way around. 

Given the exponent a_ , by series expanding the equations of motion one can determine 
the coefficients gi,hi, and all higher coefficients, in terms of /i, which is undetermined. 
However, f\ can be set to any value by rescaling the coordinates r, t, x. This reflects the 
physical fact that only ratios of dimensionful quantities are meaningful. We can therefore 
set /i to an arbitrary constant value (the sign is important however) but should make sure 
to only compute and plot dimensionless quantities. With the series expansion at hand we 
can proceed to numerically integrate to smaller values of r. A typical result is shown in 
figure [2] below. In the plot we see how the thermodynamic quantities of the fermion fluid 
flow from their constant Lifshitz values at large r to zero at the star radius r = r s . Note 
that the IR region of the spacetime, large r, has a finite volume in the radial direction. 

1.4 r 
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5 10 15 20 
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Figure 2: From bottom to top, the pressure, energy and charge density distributions for an 
electron star with z = 2 and m = 0.36 (corresponding to (3 ~ 20). The boundary of the 
star is r = r s . Recall that the boundary of spacetime is at r = while r — > 00 is the deep 
IR. In the IR the thermodynamic quantities tend to their constant Lifshitz values. 

The boundary of the star occurs when the local chemical potential is not large enough 



to populate the local Fermi sea. Thus from (2.13) 

h(r,) 
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3.3 Matching onto Reissner-Nordstrom and thermodynamics 

Outside of the electron star p = p = a = and the solution must become the (planar) 
Reissner-Nordstrom-AdS spacetime. This solution has 

2 2/o2 2 

f= C -2-Mr+ r -^-, g = ^j, h = p-rQ. (3.11) 

The four constants of integration {c, M, Q, fb}, to be related to boundary field theory quan- 
tities shortly, must be fixed by matching {f,g,h,h'} at r = r s . The perhaps unfamiliar 
constant c is necessary because the normalisation of the time coordinate has been fixed 



already by our choice of f± in (3.8) in the star interior. We could choose fx such that c = 1, 
but this will not be necessary so long as we consider dimensionless quantities. 

It is physically instructive, mimicking the standard astrophysical description of neu- 



tron stars, to let {c,M,Q} become functions of r and parametrise the solution by (3.11) 
throughout the spacetime. It is then a short exercise to show that the functions M(r) and 
Q(r) obey 

(rQ(r))' = c(r) J°° ^^a(s)ds, (3.12) 

These identities are valid for any ideal fluid and do not depend on the specific equation of 
state (although we do use the zero temperature thermodynamic relation —p = p — ah/yff). 
The above integrals show how the charge Q and energy M — \tQ 2 enclosed within a given 
radius are determined respectively by the charge density of the ideal fluid and by the sum 
of the energy density of the fluid and the energy in the electromagnetic field. 

By evaluating the previous expressions at the boundary, r = 0, we obtain formulae for 
the charge and energy densities of the dual field theory 



Q = Q(0) = c y°° ^fi-&(s)ds, (3.14) 

* - • <-> 

These quantities are densities with respect to the two boundary spatial dimensions, while 
a and p were densities with respect to the bulk three spatial dimensions. Because the UV 
theory is a relativistic conformal field theory in 2+1 dimensions, we must have that the 
pressure and energy are related by E = 2P. Furthermore, in the grand canonical ensemble 
the free energy ti, = —P. It follows from the zero temperature thermodynamic relation 

-P = E-pQ, (3.16) 
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that we must have 



E 



'-(lQ. 



(3.17) 



In these expressions fi is the chemical potential of the dual field theory, see e.g. [SJ, not to 
be confused with the local bulk chemical potential. 



The thermodynamic identity (3.17) can be checked numerically. In deriving this result 



analytically we were lead to the following useful observation: on our solutions 

2rhti - 2/ - rf = . 



(3.18) 



We can show this in two steps. Firstly we note that (3.18) is true on the Lifshitz solution 



(3.1) satisfying (3.2). Secondly, by differentiating the expression in (3.18) and using the 



general equations of motion (2.6) - (2.9) together with — p = p — ah/ ' \fj we can show that 



it remains zero along radial evolution. Now using (3.18) it is possible to integrate by parts 



in (3.15) and derive the identity (3.17). 



First integrals of the equations of motion like (3.18) are common in gravitational back- 



grounds and typically implement the isentropy of the classical gravitational flow. In our 
case we are at zero temperature and the entropy is zero. In practice, we can replace the 



second order equation (2.17) by the first order relation (3.18). Thus we have reduced the 



equations of motion to three first order equations (albeit still involving a first derivative 
squared) . 

To compare the different electron stars at different values of z and m, an instructive 
variable to consider is the dimensionless ratio of the total energy and charge (densities). 
This ratio can also be compared to the value for extremal Reissner-Nordstrom black holes 
with no fermionic hair 

^27 M 2 



1 



(extremal Reissner-Nordstrom) . 



(3.19) 



32 C Q3 

The lower this ratio, the more efficiently the solution is able to carry the charge Q. The 
ratio is shown in figure [3] for various electron stars as a function of the IR critical scaling 
exponent z for different values of the fermion mass rh. 

In figure [3] we see that the mass of an electron star at fixed charge is always lower than 
the corresponding extremal black hole with the same charge. The stars are therefore ther- 
modynamically preferred. This might appear surprising as extremal black holes are often 
the lightest charged objets in the theory with a given charge, being made of 'pure charge' 
in some sense. The situation here is quite analogous to that of holographic superconduc- 
tors [mm]. Furthermore, we should note that in an extremal black hole background, the 
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Figure 3: Dimensionless ratio of the electron star mass to its charge, normalised such that 
the ratio is unity for an extremal black hole. The three curves correspond, from top to 
bottom, to masses rh = 0.7, rh = 0.36 and rh = 0.07. 

local chemical potential h/yff — > 1 at the horizon. This can be seen by substituting the 
extremality condition ( |3.19 ) into the metric (3.11). For any fixed fermion mass rh < 1, this 



chemical potential becomes greater than the fermion mass before the horizon is reached, 
and therefore a fermion density is induced. The extremal black hole is thus never a solution 
to the equations of motion in the range rh < 1 that we are considering. Alternatively, 
one could say that the extremal black hole is a solution with an unstable vacuum for the 
fermion field, in which the fermion states with energies between m and fi\ oc . are not popu- 
lated. Consistently with previous remarks, we see that the extremal black hole behaviour 
emerges as rh — > 1 or as z — > oo. 



4 An action for charged ideal fluids 

While the equations of motion developed in previous sections are sufficient for many pur- 
poses, an action principle often plays a useful role in the holographic correspondence. In 
this section we pause in our main development to describe an action that recovers all of the 
equations of motion we have used above. 

A simple action for neutral ideal fluids coupled to gravity was formulated by Schutz [39] . 
We will start instead from an 'off-shell' form of the Schutz action, following |50J and more 
closely [51] . coupled to the Maxwell field along the lines of |52j . Off-shell refers here to the 
treatment of a constraint. 
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For simplicity we begin with the simplest case of 'non-rotating' fluids^] at zero temper- 
ature. The generalization to thermal rotating fluids is straightforward, as we will indicate 
below. An action functional describing non-rotating charged ideal fluids at zero tempera- 
ture, minimally coupled to gravity, is given by 

S = J (ftx^f^g (£ E ins. + Avlxwl. + £ fluid ) > 



where 



and 



^Eins. - 7^2 ( R + £2 ) ' £mxw1 - = "^^^ ' 



Aluid = -p{<r) + <Ju a (d a <j) + A a ) + X(u a u a + 1) . (4.3) 

We will see shortly that u a , p and a are the four velocity, the energy density and the charge 
density of the fluid respectively, A is a Lagrange multiplier, and is a 'Clebsch' potential 
variable associated with the fluid velocity. As previously, we have set the charge of the 
fermion to be unity and thereby conflated the charge and number densities. Clearly (j) must 
shift under a gauge transformation in order for the action to be gauge invariant. 



We now proceed to derive the equations of motion from (4.1 ). The variables with respect 
to which we vary the action functional are A, <r, 4>, the covariant velocity u a , and finally the 
gauge potential A a and the metric tensor g a f,. The variation with respect to the first three 
of these variables yields the following equations 

5X: \u\ = -l, (4.4) 

6a: p'(a)=u a (d a( / ) + A a ), (4.5) 
6(f>: V a (au a ) = 0. (4.6) 



Equation (4.6) is the continuity equation for the fluid current vector 

J a = au a , (4.7) 



while (4.4) is just the statement that the fluid four velocity should be timelike. On physical 



grounds we can identify the left hand side of equation (4.5) as the local chemical potential 



f,(a)=p'(a), (4.8) 



2 Because our fluid is charged, by non-rotating we will mean that e abcd ut (y c Ud + ^F c d/fi) = 0. In 
practice this condition consistently restricts the degrees of freedom of the fluid to a 'potential flow', as will 
become manifest below. 
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and introduce the fluid pressure p through the thermodynamical equation 

p(a) = -p(a) + ap(a). (4.9) 
The previous two formulae are simply useful definitions insofar as the equations of motion 



are concerned. We can rewrite (4.5) in the form 



(4.10) 



p = u a {d a (t> + A a ) 

Note that the fluid chemical potential p is gauge invariant. 

Next comes the varying of the action with respect to u b . This leads to 

Su a : a(d a <p + A a ) + 2\u a = . (4.11) 

Multiplying the previous equation by u a , the Lagrange multiplier A is determined as 

up 1 



Thereby we also obtain 



x=--f = -( P + p). 



d a <j> + A a 



(4.12) 



(4.13) 



This is the so-called 'velocity-potential representation', which is implemented directly in 
'on-shelF variational formulations. 

Varying the action with respect to the gauge potential gives the Maxwell equations 



\/ b F ab = e 2 J a , 



(4.14) 



where the current was defined in (4.7) above. Finally the Einstein equations for the geom- 
etry, 



1 



1 



Rab ~ ^9abR ~ j^gab = ^ \F ac F b c - ^g ab F cd F cd ) + k 2 T^ id , (1.1 5 ) 



are obtained upon varying the action (4.1) with respect to g a b- The energy-momentum 
tensor of the fluid is given by 
2 5 



rpab 
1 fluid 



-g 5g a b 



J d 4 xy/=g £ fluid = <f % uid - 2u ( ~ a [a(d b U + A b ^) + An")] . (4.16) 



Using the previous equations of motion, the on-shell energy-momentum is found to take the 
standard form for an ideal fluid 



T£ id = (p + p)u a u b +pg a \ 



a-„.b 



..ob 



(4.17) 



We should check that we have indeed recovered all the equations of the previous section 



from the full action (4.1). The Einstein-Maxwell-fluid equations (2.1) - (2.3) have been 
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obtained explicitly. The equation of state p{a) was denned implicitly through equation 



(2.10), by elimination of \i. We can easily check that (2.10) furthermore implies fi = p'{<y) 



as required. The ansatz we made in section [2] for the metric, Maxwell field and fluid velocity 
corresponds to setting (f> = in (4.13) with the local /i consequently given by ( |2.12| )plThus 
we see that all the equations used in previous sections correspond to a solution of the 



equations of motion following from the action (4.1). 

It is now easy to obtain the on-shell form of £ nu id- This is required for instance to 
evaluate the free energy of the electron star. Substituting the equations of motion into the 



action (4.3) gives 



£Z7* hen =P- (4-18) 

Thus the on-shell Lagrangian for the ideal fluid is simply its pressure. Evaluated on our 



electron star ansatz (2.4), one can then verify using the equations of motion (2.6) - (2.9) 



(4.19) 



that the full Lagrangian in (4.1) is a total derivative 

£ on- S hell _ L 2 d f - 2hh> ^ 

k 2 dr 2r 2 \pfg 

The free energy is given by the Euclidean action evaluated on shell. We have just seen that 
the bulk action becomes a boundary term. In order to obtain a finite answer we must add 
the boundary Gibbons-Hawking term and intrinsic counterterms. We will not describe this 
standard process in detail, see e.g. [8J. The solution near the boundary takes the form 



(3.11), this is all that is needed to evaluate the action on shell as there is no contribution 



from the IR Lifshitz endpoint of the integral. The upshot is that the free energy density is 



O = M - jXQ . 



(4.20) 



as we assumed in (3.16) above. 



We can also use the action (4.1) as a starting point for Lifshitz holography, with the 



UV given by (3.1) rather than AdS±. Lifshitz holography requires additional boundary 



counterterms. These are most conveniently packaged [531 154j as a series in powers of \d(j) + 
A\ 2j t- e ^ , as this combination vanishes on the Lifshitz background and it becomes apparent 
that only a finite number of such terms are necessary. 

3 As usual in holographic setups, we fixed a gauge ambiguity by requiring At to vanish at the horizon. 



Adding a nonzero constant d t (f> in the formula for the gauge invariant local chemical potential (4.10l would 
result in a non-Lifshitz invariant IR. A radially dependent <f> would take us outside of our ansatz, introducing 
radial fluid flow. This choice of requiring fj, to tend to a constant in the IR can presumably be thought of as 
a choice of fermion vacuum that is regular at the IR 'horizon'. Thanks to Tom Hartman for discussions of 
this point. 
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If, following Schutz |49j . we wished to obtain the correct equations of motion for the 



fluid as well as the geometry from the 'on-shelP Lagrangian (4.18), we would firstly need to 



bring the equation of state of the fluid p = p{a) into the form p = p(p). This is achieved 



through (4.9) which may be viewed as the usual Legendre transformation. Subsequently we 



can take the norm of (4.13) to express the pressure in terms of the Clebsch potential and 



the gauge potential A a . This leads to the final Schutz form of the fluid Lagrangian 

= P(») = P(\d<P + A\) . 



(4.21) 



This action should be varied with respect to <j), A a and g a b- This form of the action shows 
most explicitly that coupling a charged ideal fluid to gravity without rotation or temperature 
is equivalent to coupling to a Stiickelberg field. 



Finally, we should sketch the straightforward generalization of the off-shell action (4.1) 



to describe charged rotating and finite temperature ideal fluids coupled to gravity. More 
details can be found in the papers we referred to above. One first introduces two pairs of 
new potential- variables, (a, (3) and (s,9). The first pair will account for the fluid rotation. 
In the second pair, s will become the fluid entropy density, while the variable 9, the so-called 



'thermasy', will be responsible for the fluid temperature. Accordingly, in the action (4.1), 
the fluid Lagrangian is promoted to the following 



£fluid = ~p((T, s) + au a (d a cp + A a + 0d a s + ad a [3) + \{u a u a + 1) . 



(4.22) 



The equation of state has been enhanced to include an entropy dependence, p(a,s). The 



equations (4.4) and (4.6) will remain the same, while equation (4.11) is replaced with 



a{d a + A a + 9d a s + ad a p) + 2\u a = . 



(4.23) 



There is a similar modification to equation (4.10). Further, we have the following equations 



from varying the action with respect to the new potential variables 



u a d a s = ; u a d a f3 = ; u a d a a = ; u a d a 9 = -T 



1 dp 

ds 



a 



(4.24) 



in which T denotes the fluid temperature. Similarly to above, the resulting Schutz form of 
the action is found to be 

p(fi,s), (4.25) 



p Schutz 
-'-fluid 



where now the chemical potential in terms of the potential-variables and the gauge field is 



p = \d(j) + a df3 + 9 ds + A\ . 



(4.26) 
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5 Electrical conductivity 

To compute the frequency dependent electrical conductivity at zero momentum, a(ui), we 
need to perturb the backgrounds of the previous sections. We clearly need to perturb the 
vector potential A x , as this is the field dual to the electric current. At zero momentum, 
these perturbations source perturbations of the metric component gt x and the velocity u x . 
Specifically, if we take the perturbations to have time dependence e~ lujt , so that 

A x = —5A x (r)e-^\ g tx = L 2 5g tx (r)e~ iut , u x = L 5u x (r)e~ iult , (5.1) 

then the linearised Einstein-Maxwell equations about the above backgrounds are solved if 
the following three equations are satisfied 

aSA x + (p + p)5u x = 0, (5.2) 

Sg' tx + -5g tx + 2h'5A x = 0, (5.3) 
r 

5 ^ + \ij~ ^) 5A '* + j( 6 9L + ls9t^+gv5 Ux + u; 2 ^5A x = 0. (5.4) 

It is immediately clear that we can eliminate 5gt x and Su x from the above equations to 
obtain a single equation for 5A X 

SA X + - 2 {j--) SA X + {—-— r —) SA9 = 0. (5.5) 



The structure of this equation is similar to that arising in holographic superconductors 
There is a 'mass' term due to screening in this case rather than electromagnetic symmetry 



breaking. The rightmost term in the equation (5.5) is due to the coupling between metric 
and Maxwell fluctuations. It is also a mass-like term and leads to an infinite DC conductivity 
because the medium has a net charge density and is translation invariant |17} IT8] . 



The equation for the perturbation of the fluid velocity (5.2) is just 5A X + [ibu x = 0, 



which is again compatible with the irrotational form (4.13) with Sfj, = 6(f) = 0. It is possible 



that modes with finite momentum k will excite the scalar degree of freedom of the fluid. 
5.1 The conductivity at low frequencies 



In the Lifshitz IR background (3.1) we can solve the equation (5.5) for the 5A X fluctuations 



analytically in terms of a Hankel function 



6AW) = ^H$(&^ . (5.6) 

In deriving this formula we imposed ingoing boundary conditions at the Lifshitz 'horizon' 
(i.e. 5A X ~ e+^ffoo r z /z as r _^ We also used the algebraic equations (3.2) giving 
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Socn^oo in terms of z,/3,rh in order to simplify the index of the Hankel function. Similar 
Hankel functions were found for fluctuations of a vector field in a Lifshitz background in e.g. 
\45\ [37] 114] . In fact, the use of Hankel functions here is overkill, as the previous expression 
can equivalently be written as an oscillating exponential multiplying a polynomial 

SAW + (5 .7) 



We have changed the (unobservable) overall normalisation relative to (5.6). 



The solution for the perturbation in the Lifshitz region (5.7) will hold for the full solution 
in the 'near' region defined by r/i > 1. At zero temperature £i is the only energy scale in 
the problem. This condition simply means that we can use the leading order metric near 
the horizon. At low frequencies, uj -C fi, the near region has an overlap with the 'far' region 
defined by wr 2 /^ -1 <C 1. In the far region, away from the non-analytic ingoing boundary 
conditions, we can set co = in the Maxwell equation in order to compute to leading order at 



low frequencies. We will now proceed to match the near solution (5.7) to the far region and 
obtain the conductivity to leading order at low frequencies. This computation is essentially 
identical to that appearing in e.g. [S3 [33 SS] ■ 

The conductivity of the dual field theory is computed at the AdS& boundary. Near the 
AdSi boundary, r — > 0, the Maxwell field perturbation behaves as 

SA X = <5A(°) + rSA^ + ■■■ . (5.8) 

The 'dimensionless' conductivity is then given by (e.g. [8]) 

- — 2 * C ^ /r n\ 

a = e a = M . (5.9) 

u SA X °> 

The extra factor of c compared to more common expressions is due to the normalisation of 



the metric in 3.11 We use the same symbol for the conductivity of the dual field theory 
and the charge density in the bulk. Hopefully the context will make it obvious which we 
are referring to. 

The following flux is independent of the radial position r 



T = i^/Th {SA X 5A' X - 5A X 5A' X ) . (5.10) 



Constancy of this quantity follows directly from the equation (5.5). Evaluating near the 



boundary and using (5.9) gives 

T = 2\5A^\ 2 ooRe&(uj) . (5.11) 
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Evaluating using the near horizon solution (5.7) gives 



F = 2uj. (5.12) 
Equating the previous two formulae we see that to obtain the real part of the conductivity 



it only remains to compute the 5A X ' obtained by matching onto (5.7). Expanding the near 



solution (5.7) into the matching region uir z fi z <S 1 gives 



5A X = ——. (5.13) 

i/oo 

As we noted above, in the far region we can drop the to dependence in the Maxwell equation 



(5.5) because ui <C jl. It follows that the oj dependence in the previous formula will remain 
the same all the way out to the boundary, leading (generically) to 5A^ oc cj _1 . Putting all 
of these facts together then leads to the conclusion that to leading order at low (but finite) 
frequencies 

Reo-(w) oc oj 2 . (5.14) 

This is precisely the same low frequency behaviour for the conductivity in a dual geometry 
with an IR Lifshitz region as that obtained in |37| 05] . The same behaviour also occurs 
at extremal black hole horizons [56 1 . The physics behind the emergence of the Lifshitz 
region is apparently distinct in these various cases; this seeming universality in the electrical 
conductivity remains to be properly understood. 



As well as the real part (5.14), we can anticipate that the imaginary part of the con- 
ductivity will have a pole as u — > corresponding to a delta function at u = in the 
real part. Thus in fact Re<r(a;) oc 5(uj) + oj 2 . This delta function is due to the fact that 
the system is translationally invariant and carries a net charge. When excited by a time 
independent electric field, the whole system is accelerated, leading to a current that cannot 
be relaxed [171 [T8] . We will see the pole in the imaginary part of the conductivity shortly 
in our numerics. 

Before moving on to compute the full conductivity numerically, we should make a remark 
about the result that Rea(w) oc 5(oj) + oj 2 . Namely, that there are two reasons why the 
to 2 conductivity (which will presumably translate into a T 2 low temperature dependence of 
the DC conductivity as in |57j ) should not be taken overly seriously as an 'experimental' 
feature of this model. Firstly, in practice, disorder or other physics will smear out the 
delta function into a Drude peak. The magnitude of this peak could easily dominate, e.g. 
the temperature dependence of the DC conductivity. Secondly, putting aside the Drude 
peak, an oj 2 or T 2 dependence of the conductivity is very weak and corresponds to a huge 
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resistivity, going like w -2 or T~ 2 . It will be very easy for any other conduction channel 
to short-circuit this classical contribution. For instance, one loop processes in the bulk 
involving charged fermions, analogous to those of |58j, will likely have a resistivity that goes 
to zero at low temperatures or frequencies. The semiclassical expansion will therefore break 
down at sufficiently low frequencies at which non-classical conduction becomes favoured. It 
is clearly of interest to investigate the one loop physics of our electron star backgrounds. 



5.2 The full conductivity 

In this subsection we compute the full frequency dependent conductivity. To do this we 



must numerically solve the differential equation (5.5 ). We can integrate out from the horizon 



to the boundary and then read off the conductivity using (5.9). The ingoing solution at the 
horizon, r — > oo, on a general background takes the form 

1/2/ 



6A* 



Ju>g 1 l 2 (r*/z + r z + a -(g 1 -f 1 )/2(z+a-))( 1 + ^z+Icl- + . \ 



(5.15) 



Here a_ < is the exponent of the IR irrelevant mode in (3.9) while g\ and f\ are the 



coefficients appearing in (3.8). These are determined, along with # and higher order terms, 



by series expanding the equations of motion. From the definition of a_ in (3.9) and the 



lower bound on z in (3.7) we can show that z + a_ > and z + 2a_ < 0. This is why the 



first power must be kept in the exponent while the second can be expanded as r — > oo. 



Starting from the series expansion (5.15), we can numerically integrate the Maxwell 



equation (5.5) out to the electron star boundary r s . Outside of the electron star, we need 



to solve the equation in the Reissner-Nordstom-AdS background (3.11). This must also be 



done numerically, with the value and derivative of the fluctuation 6A X matched across the 



electron star boundary. Integrating out to the boundary then gives the conductivity (5.9) 



We performed these integrations using NDSolve in Mathematica. The resulting real and 
imaginary parts of the electrical conductivity as a function of frequency are shown in figure 
H below. 

In the plots of figure [I] we see the expected soft gap a ~ uj 2 of the (dissipative) real part 



of the conductivity at low frequencies anticipated in (5.14) above. The divergence of the 



imaginary part indicates the presence of a delta function in the real part, via for instance 
the Kramers-Kronig relations. As noted above the divergence is due to the combination 
of a net charge density and translation invariance (i.e. no impurities or lattice). At large 
frequencies the real part tends to a constant. This follows from the fact that conductivity is 
dimensionless in 2+1 dimensions and that the UV completion of our boundary field theory 
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Figure 4: The zero temperature real and imaginary parts of the electrical conductivity as a 
function of frequency. From left to right in each plot: {z = 3, to = 0.7}, {z = 2, to = 0.36} 
and {z = 1.5, to = 0.15}. The real part also contains a delta function at oo = 0. 

is a conformal fixed point with no inherent scale. 

The z = 3 plot of figure [4] is already quite similar to the zero temperature limit of the 
z = oo Reissner-Nordstrom case, which is a charged black hole rather than an electron star. 
See e.g. [8]. The only qualitative effect of the electron star (i.e. lower mass fermions and 
lower IR scaling z) seems to be to smoothen out the transition from the IR uj 2 scaling to 
the constant high frequency behaviour. The coefficient of the pole in the imaginary 

part of the conductivity is proportional to Q/fi 2 oc Q 3 /E 2 , again see e.g. [Sj. It is therefore 
consistent with figure [3] that we see that the pole is stronger at lower fermion mass. At 
lower fermion mass, the electron star has a larger charge at fixed chemical potential. 

6 Final comments 

One important objective of holographic approaches to condensed matter is to characterise 
possible (computationally controlled) IR fixed point behaviour that falls outside of the 
Landau Fermi liquid paradigm. As with the Fermi liquid itself, this is in the first instance 
a question about universal low energy physics, not about the UV physics (be it electrons in 
a lattice or some cousin of J\f = 8 super Yang-Mills theory) . This perspective suggests that 
in the bulk one should focus on the near horizon region of the geometry, as argued most 
explicitly in \59\ 160] . However, the role of a finite charge density is subtle in this regard]^] 
The charge density itself, or the chemical potential, is a UV quantity that is specified at 
4 We would like to acknowledge helpful discussions with Diego Hofman on this topic. 
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the boundary of the bulk geometry. The deep IR Lifshitz solution, in our case for instance, 
does not immediately 'know' what the value of this charge density is. In fact, the electric 
field is zero at the Lifshitz 'horizon' and grows as one moves out towards the boundary of 
the electron star. Therefore, if we wish to understand how the emergent IR Lifshitz scaling 
is related to the fact that we are considering a system at finite density, we need to connect 
statements about the UV and IR physics. In a Fermi liquid, the connection between UV 
and IR physics is achieved via the Luttinger theorem |61| [62] . This theorem states that the 
volume enclosed by the Fermi surface is determined by the average particle number. When 
the low energy effective field theory is Fermi liquid theory, the theorem essentially reduces 
to counting charged states in the UV and IR [63J. A pressing open question in applications 
of holography to condensed matter systems is to formulate a useful holographic analogue 
of this theorem. 

Partially motivated to obtain an arena where an interesting low energy scaling geometry 
could be seen to emerge from a finite density system, in this paper we have constructed 
zero temperature 'electron star' solutions in asymptotically AdS± spacetime, building on 
results in [T3J. These geometries are solutions to the Einstein-Maxwell-charged ideal fluid 
equations of motion. The structure of the solutions is a 'domain wall' flow from a 'near 
horizon' Lifshitz geometry to AdSi at high energies. The essential physics of the flow is that 
the Maxwell field becomes screened in the Lifshitz region by the charged fluid. In several 
regards our solutions are qualitatively (and quantitively!) similar to the zero temperature 
holographic superconductors of [361 GZ] an d the extremal dilatonic black holes of [IS] . 

Another motivation of our work was to obtain gravitational duals in which the full 
charge density was manifestly 'fermionic'. The electron star is literally a Fermi surface 
that is inhomogeneous in the bulk radial direction. This may or may not be an important 
ingredient in formulating a holographic Luttinger-like theorem. There are various immediate 
questions to be explored in this regard. One should compute the momentum dependence of 
the conductivity to look for Fermi surface related signatures. Also, upon adding a magnetic 
field to the system the electron star should show quantum oscillations already at a classical 
level, unlike charged black holes for which quantum oscillations are only present at one loop 
order and reveal a 'small' Fermi surface \25\ 164] . 

There are various directions in which our work could be extended at the level of general- 
ising the solutions we have presented and studying their physics. Upon adding interactions 
the stars will likely have instabilities such as Cooper pairing instabilities along the lines of 
|27] , One should explore the effects of changing the equation of state. There may be cir- 
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cumstances in which clumping instabilities arise. It should be straightforward to place the 
system at finite temperature. It is possible that a phase transition to a black hole solution 
occurs at some finite temperature, analogously to the transition in the neutron stars of [43] . 
A more challenging question is to move away from the ideal fluid limit. One would like 
to solve the Dirac equation in an unspecified background, populate the low lying states up 
to some chemical potential and then self-consistently solve the Einstein-Maxwell equations 
together with this quantum source. 
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